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The sensitivity of polarization observables in the reaction NN — > NN(p 
slightly above the threshold is studied with respect to the details of the one- 
boson exchange model and a possible admixture of hidden strangeness in the 
nucleon. It is shown that the finite-energy predictions differ strongly from the 
threshold predictions. A measurement of the beam - target asymmetry and 
the <j) decay anisotropy can help to disentangle the role of various reaction 
mechanisms. 
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I. INTRODUCTION 

The investigation of the reaction NN —>■ NN(p is interesting for several reasons. First, 
the elementary total cross section is an important input for the calculation of the 
production in heavy-ion collisions ||. In this case one might expect some change of the 
width || due to the coupling to the decay channel — > K + K~ and peculiarities according 
to the in-medium modification of the kaon properties Indeed, such measurements are 

under way with the An detector FOPI at SIS in GSI/Darmstadt 0. The electromagnetic 
decay channel — > e + e~ will be investigated with the spectrometer HADES also in 
GSI. A threshold-near measurement of the total cross section of the reaction pp — > pp<f) has 
been performed at SATURNE || and precision measurements of polarization observables 
are envisaged with the ANKE spectrometer at the cooler synchrotron COSY in Jiilich ||. 

Second, the meson is thought to consist mainly of strange quarks, i.e. ss, with a 
rather small contribution of the light u, d quarks. According to the OZI rule (cf. |T(J for 
references) the production should be suppressed if the entrance channel does not possess 
a considerable admixture of hidden or open strangeness. Indeed, the recent experiments on 
the proton annihilation at rest (cf. |TIJ for a compilation of data) point to a large apparent 



violation of the OZI rule, which is interpreted Jl0| as a hint to an intrinsic ss component in 



the proton. However, the data can be explained by modified meson exchange models [11 



without any strangeness content of nucleon or OZI rule violation. On the other hand the 



analysis of the 7tN sigma term [12] suggest that the proton might contain a strange quark 



admixture as large as 20%. Thus this issue remains rather controversial. 



Therefore it is tempting to look for other observables ||T0| , p!3| that are directly related to 



the strangeness content of the nucleon. In Ref. |T3[] it is shown that polarization observables 



of the photoproduction are sensitive even to a small strangeness admixture in the proton. 
Other investigations concern the polarization observables in pp — > pp<p reactions at the 
threshold fltl| , |15|| . It is found Jl5] that spin and parity conservation arguments result in a 



precise prediction for the beam - target asymmetry 
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^ BT da(Si = 1) - da (Si = 0) 



dcr(5i = 1) + da(Si = 0) 1 J 

for the pp — ► pp0 reaction at the threshold, C^, T = 1, where is the total spin in the entrance 
channel. It is further claimed [15| that the elements of the decay density matrix p rr i (for 
notation see below) amount p 00 = 0, p±\±i = |. Real experiments, however, are performed 
at some finite energy above the threshold, and therefore the problem arises whether the 
threshold predictions are modified. More generally one should investigate the sensitivity 
of polarization observables to the strangeness admixture in the proton. An analysis of the 
threshold-near behavior of the unpolarized cross section is a necessary prerequisite, of course. 

In the present work we employ a one-boson exchange (OBE) model to consider these 
mentioned issues. In Section II we present expressions for the kinematics, cross sections and 
spin density matrix. In Section III we describe our OBE approach. Section IV is devoted to 
parameter fixing. The threshold limit is analyzed in Section V. The results at finite energy 
above threshold are discussed in Section VI and summarized in Section VII. 



II. KINEMATICS AND OBSERVABLES 

The corresponding cross section of <fr production in the reaction a + b — > c + d + <p, where 
a, b and c, d label the incoming and outgoing nucleons (cf. Fig. 1), is related to the invariant 
amplitude T as 

*>= _. .. - 1 PV^^W-P/)- (2) 



2(2tt) 5 ^/s(s - 2E c 2E d 2E 4 



where p n = (E n , p n ) with n = a,b,c,d and q = (E^, q) are the four-momenta of the nucleons 
and the meson in the center of mass system (c.m.s.), respectively. Hereafter 9 denotes the 
polar cf) meson angle, s = E a + Eb is the total c.m.s. energy, and Pij is the total 4-momentum 
of the initial or final states. We use a coordinate system with z || p a , y || p a xq. Among the 
five independent variables for describing the final state we choose E^, and Q c . Then the 
energy E c of particle c is defined by 



AB±C^B*-M%(A*-C*) 
c = A 2 -C 2 ' ^ ' 

with A = 2(t/s — E ( / > ), B = s — 2E ( f ) T/s + Ml, C = 2qcos9 qPc , where 9 qPc is the angle between 
q and p c determined by cos9 qPc = cos 9 cos 9 C + sin 9 sin 9 C cos ((p — (p c ); the unindexed polar 
and azimuthal angles 9 and <fi belong to the vector q. In calculations one has to use both 
solutions of Eq. ([3D according to the considered kinematics. Finally, the fivefold differential 
cross section reads 

d b o 1 i T|2 qvl m 



dE^dn^dn, 8 ( 2vr ) s ^ s ( s _ 4M 2 ) | a Pc + ce, 

The total cross section is found by integration of over the available phase space 
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One can consider the asymmetry Eq. (1) with integrated cross section or differential cross 
sections. In the general case the asymmetry depends also on E^, and Q c . 

We consider the decay angular distribution W of the <f> meson in its helicity system where 
the <fi is at rest. The quantization axis z' is chosen along the meson velocity in the overall 
c.m.s., and the y' axis is parallel to the vector p a x q. The decay angles 6, <£> are defined 
as polar and azimuthal angles of the direction of flight of one of the decay particles in the 
meson's rest frame. The decay angular distribution then reads 

dN 

W(cos 9, $) = M ( r ; e > $ ) Prr> M*(r'] O, $), (5) 



dcosOd® 



where M(r; 6, $) is the decay amplitude, M(r; 6, $) = M d x r * Xa _ x ($, 6, -$), with r as 



the helicity of meson; A Qi/ 3 denotes the helicity of decay particles and |M | 2 is proportional 
to the decay width; D\ , stands for Wigner's rotation function. The angular distribution is 
normalized as 



W(cos0,$)dft = 1. (6) 
The density matrix p rr t is expressed in terms of the production amplitudes 

(7) 



Ernr T*r 



2^r,/3 1 (3 1 p 

where j3 denotes a set of unobserved quantum numbers. 

Assuming that in — > K + K~ and — > e + e~ decays \ a — \p = (i.e. A Q — A^ = ±1), 
and using the explicit form of the rotation function we get 

W K+K ~ (cos 9, $) = A [ i( Pll + p_!_!) sin 2 + poo cos 2 

47T 2 

H — W— Repio + Rep_i ) sin 20 cos $ H — Wlmpio + Imp_i ) sin 20 sin $ 
v2 v2 

- Repi_i sin 2 cos 2$ + Impi_i sin 2 sin 2$ ] , (8) 
W e+e ~ (cos O, $) = J- [ h Pu + p.^l + cos 2 0) + poo sin 2 

57T 2 

H ^(Repio — Rep_i ) sin 20 cos $ Wlmpio + Imp_io) sin 20 sin $ 

+ ^Repi^i sin 2 cos 2$ - ^Impi_i sin 2 sin 2$ ] . (9) 

The decay distributions integrated over the azimuthal angle depend only on the diagonal 
matrix elements, 

W^^cosO) = |(1 - Poo) (l + B K+R - cos 2 0) , 
W e+e ~ (cos 0) = |(1 + poo) (l + B e+e ~ cos 2 e) , (10) 



where we use the normalization condition p n + p_i_i + poo = 
anisotropies B 



1 and introduce the decay 



^ + ^ = -(l-3p o)/(l-Poo), 
5 e+e ~ = (l-3p 00 )/(l + Poo). (11) 



III. ONE-BOSON EXCHANGE MODEL 



Similar to previous analyses of the bremsstrahlung in pp reactions |16| we use here the 



OBE model to parameterize the interaction in the reaction iViV — > NN<f> by the set of 
relevant tree level diagrams with one meson exchange line. These reaction channels are 
depicted in Fig. 1. In diagram la the is produced in an internal conversion in a 0p7r vertex, 
while diagram lb is the contribution of the "bremsstrahlung" due to a finite coupling of 
the to the nucleon current, which might arise microscopically from the coupling of the 
to the kaon cloud around the nucleon. Diagram lc depicts a more exotic channel, namely 
the shake-off reaction by the internal virtual meson flow analog to the photoproduction 



bl] , |TT|| . It represents the emission from the vertex. (As alternative one could consider the 
ss knock-out, but we assume here that the corresponding meson-ss-0 coupling constants are 
small.) Within the OBE model the diagrams la - c cover basically all reaction mechanisms. 
One expects that diagram la represents the dominating contribution to the OBE mechanism, 
and we call this process hereafter the conventional or pure OBE contribution. 

A. Emission from internal meson line 

This part of our OBE model is similar to the previously employed model ||. The meson- 
nucleon and the 0p7r interaction Lagrangians read in obvious standard notation 



C M nn = -ig^NNNysTTrN - g pNN {N^rNp" - la* v rNd uPll J , (12) 

Cfa = 9^ ^ d^ u Tr(d aPf3 n), (13) 

where Tr(p7r) = p°7r° + p + 7r~ + p _ vr + , and bold face letters denote isovectors. All coupling 
constants are dressed by monopole form factors Fi = (A? — mf)/(Aj — kf), where hi is the 
four-momentum of the exchanged meson. 

The total invariant amplitude is the sum of 4 amplitudes 

T r a = E {£T r [ab; cd] + £T r [ab; dc] + £T r [ba; dc] + &T r [ba; cd}) (14) 

wi th £p P = £p P = -tip = -£p P = 1> £pn = £pn = -1> £pn = £pn = 2 - Th e last two terms 
stem from the antisymmetrization or from charged meson exchange for pp or pn reactions, 
respectively. The first term in Eq. (|i"4|) for the pp reaction reads 

T^n cmsma m b [ob\ cd] = K(l£, Tlm d m b {Pb, Vd) C c m„ (Pa, Pc) (15) 
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u(p a ,m a 



(16) 



Here e r,/3 denotes the polarization four-vector, and m a ...m d are the nucleon spin projections 
on the quantization axis, and k p = p c — p a , k n = p b — p d ; u stands for a Dirac bispinor. 
Detailed expressions of the functions n, W and £ are displayed in the Appendix. 

The amplitude in the helicity frame is obtained from the amplitude in the c.m.s. by a 
standard Wigner rotation. 



B. Emission from external legs 

The amplitude of direct <fi meson emission from the nucleon legs according to Fig. lb 
is calculated similar to the real or virtual photon bremsstrahlung [18|.[19| in the few GeV 
region. The (pNN vertex has the same structure as the pNN interaction (cf. second part of 
Eq. fljjD), i.e., it possesses vector and tensor parts. Following Ref. [[2(J we use q^nn = —0.24 
and = 0.2. 

The internal zig-zag lines in Fig. lb correspond to a suitably parameterized AW inter- 
action in the form of the Born term of the OBE model with effective coupling constants 
and cut-off parameters and may be interpreted as effective Tr,u,p,a meson exchanges. We 
would like to stress that this is an effective dynamical model which is appropriate in the 
few GeV region and which is different from the OBE model in the conventional sense. This 



model has been applied successfully to different reactions |p!8| , [19| , [16| and this encourages us 
to employ it for the (f> production too. 

The total amplitude for the process lb consists of 4 • 8 contributions and has a similar 
structure as Eq. flT3D, where T[ab; cd] now reads 

T r [ab; cd] = - g<f)NN ]T (u(p c ) [T m fij + n r 2 T m ] u{p a )) [-i D m ] (u(p d )T m u(p b )) , 

m=7T,(T,p,UJ 

l ~p 2 -M 2 N e P [lmU + l 2M N a quU M 2M N qA p 2 -M 2 N {U) 

Here we use the notation p\ = p a — q, P2 = Pc + Q, and T m and D m are effective coupling 
vertices and propagators of the two-body T matrix, respectively, 



k 2 — m 2 i' Mi/ k 2 — m 2 



T n = -ig^NNls, ^ = g aNN , = -g p/uNN ( 7a , + i^u^K ) , (18) 
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where k is the momentum of virtual meson and g m NN is the corresponding coupling constant 
which is dressed by the above form factors. The explicit expression of the total amplitude 
is rather cumbersome and will not displayed it here. The most economic way to evaluate 
this amplitude is a direct numerical calculation of the Feynman diagrams including the 
spinor algebra. This method has been checked in a number of examples (including the OBE 
amplitude in Fig. la). 



C. Emission from vertex 



The main ingredient of the <fi shake-off mechanism (cf. Fig. lc) is the assumption that 
the constituent quark wave function of the proton contains, in addition to the usual uud 



component, a configuration with an explicit ss contribution |T7|j2T| . A simple realization of 
this picture is the following wave function in the Fock space 



m 



\p) = A\[uud} 1 ' 2 ) + B{a \[[uud] 1 ' 2 ® M°] 1/2 > + a^uud} 1 ' 2 <g> [ss] 1 ] 1 / 2 )}, 



(19) 



where B 2 denotes the strangeness admixture (|A| 2 + |£>| 2 = 1), and a 2 , x = | are the fractions 
of the ss pair with spin and 1. The superscripts represent the spin of each cluster, and 
® indicates the vector addition of spins of the uud and ss clusters and their relative orbital 
angular momentum (£ = 1). Details of the wave functions in the relativistic harmonic oscil- 
lator model can be found in [FJ]. Some estimates regarding the strangeness content in the 
nucleon advocate a few percent |10|JT3|| . Otherwise, calculations based on the inclusion of 
extra kaonic degrees of freedom within a perturbative treatment give more than an order of 
magnitude smaller value of hidden strangeness fl4]j2"2|] which is, however, likely to be domi- 
nated by nonperturbative effects [ 23| . That is the reason why we consider <p bremsstrahlung 
and shake-off independently. 

We assume here as above that the exchanged effective mesons interact with the uud 
cluster as with a structureless particle and parameterize this interaction within the effective 
two-body T matrix. The ss component is considered as spectator, that means only the 
configuration with spin S sS = 1 is realized. The corresponding T matrix element is calculated 
in the rest frame of the splitting nucleon (p a = (M N , 0)) and it is expressed by 



T 



m c m d ;m a in b 



[ab : cd}^ 



E E 

=n,p,uj,a m x 



M rN^ T (m)NN^NN ( / x 

m x :m a m c m d \m x m, h KFci ydi Fa) fb ) 



(20) 



with T( m ) 



NN^NN 



as effective two-body T matrix 



4£T(Pc,Mift) = (u(p c ,m c )T m u(p' a m x )) [-zD m ] (u(p d ,m d )T m u(p b ,m b )) (21) 
and with the vertices and propagators as in Eq. ([18]). The transition amplitude reads 



M. 



m x ;m a 




jm. 



3rrij-m x 



1 

—m, 

2 c 



with 
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V(a L ) 

C m = A*Ba l ^ 1 f, (22) 

and p' a = zp a , where z ~ 0.5 is the momentum fraction carried by the uud cluster, q L 
denotes the (f) momentum in the laboratory system, is its projected unit vector in the 
circular basis, and V(\t) stands for the wave function of the relative motion normalized as 
/ <ikV 2 (k)/((27r) 3 24/(k 2 + M?)) = 1; 7<^ c is the corresponding Lorentz factor which reflects 



the Lorentz contraction in the relativistic constituent model. In our calculations we use a 
Gaussian distribution V(x) = Nx exp(-x 2 /2fi) with = 2.41 fm" 1 jl4lJl7| . 



By making use of the numerical values of the Clebsch-Gordan coefficients we get the final 
expression for the amplitude M N ^ N< ^ in the form 

V 1 m =O.±l V 2 

(777, * \ 
-j=-m a \) )d 1 m4>r (&). (23) 

The two body scattering amplitude T^J*?^^ N for each exchange meson is found straight- 
forwardly by using the functions II and £ in the Appendix. The final amplitude contains a 
sum over all exchanged mesons and consists of 2 direct and 2 exchange amplitudes for the 
pp reaction, and 2 direct and (lj, a) or 2 (n, p) exchange amplitudes for the pn reaction to 
be taken with their proper isospin factors £^,( pn ) as in Eq. flT4]). 



IV. FIXING PARAMETERS 

The parameters of the two-body T matrix for <p bremsstrahlung (Fig. lb) and shake-off 
(Fig. lc) are taken from Refs. [[I6i[l8| where quite reasonable agreement with data of different 
elastic and inelastic NN reactions is found. 

The coupling constant is determined by the 4> — > pir decay. The recent value 
r(0 — > pit) = 0.69 MeV results in g^ pn = 1.10 GeV -1 . The remaining parameters of the OBE 
amplitude for the process la are taken from the Bonn model as listed in Table B.l (Model II) 
of Ref. P4]| . The cut-off parameter A^ p7r = 2.9 GeV is adjusted by a comparison of additional 



calculations including the channels discussed above and data |25| for the n~p — > rxj) reaction. 
Note, that on this stage the result is not sensitive to the specific OBE model of the irp 
reaction, where A^ p7r is used as fit parameter, that is, one can fit the data also with any 
other OBE model by using A^ p7r as a free parameter. Turning to the pp reaction, however, we 
find that the cut-off A^ p7r represents the main uncertainty of the model. In would be perfect 
to find it from the unpolarized cross section in this energy region and after that to study 
the polarization observables. Lacking precise data one has to use some reasonable guess 
which make our finite energy prediction rather qualitative estimates than exact predictions. 
However, as we show in the next section, some of the threshold prediction are "parameter 
independent" . In any case, using our model one can easily refine our finite energy prediction 
when the needed data on unpolarized <fi production are available and fix Aj p7r . The value 
A^ p7r = 0.77 GeV is provided by an analysis of the photoproduction of vector mesons p6| , 
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and we use here this value having in mind, however, that the symmetry of the off-shell 
mesons in the (fipn vertex would suggest = A^ pn which can be employed as another 
reference value. 

Note that the OBE model describes the data |25| in a limited region of A^-p^^ = 
\fs — Mn — Mff, < 0.15 GeV. At higher energies the OBE model overestimates the data. 
(In this region one could use an energy depending suppression factor at the vertices as in 
Ref. 1). 

Another type of uncertainty is related to the total phase (rj = exp(— z<5)) of the two-body 
T matrix. In the effective theory of Refs. [18, TBI one fixes only the relative phase between 
different exchange amplitudes and, therefore, there are different possible solutions. The 
most natural one is the assumption that each exchange amplitude has the same phase as the 
pure amplitude of the process la, which is known. In this case the two-body T matrix must 
be real (5 = 0). In that instant we expect a strong interference between the radiation from 
the internal line and the external legs, because both amplitudes are almost real (at least 
in coplanar geometry), and a small interference with the radiation from the vertex, where 
the amplitude is almost imaginary. The opposite limit is based on the assumption that 
the unitary condition SS + = 1 and the relation S = 1 — iT result in a forward scattering 
amplitude which is almost imaginary and negative. The assumption that the effective T 
matrix is responsible for the elastic forward scattering constrains the common phase, at 
least at small momentum transfer (5 = vr/2). In this case one gets a strong interference 
between the OBE amplitude la and the radiation from vertex lc, and a small one between 
these and the radiation from the legs lb. The latter assumption is more reliable for high 
energy scattering, where the diffractive scattering (i.e. Pomeron exchange) is dominant; in 
near-threshold production region it has rather a methodical interest. Nevertheless, for 
completeness we will consider both cases. 



V. THRESHOLD LIMIT 

Let us now consider in some detail the threshold limit, where one can neglect terms 
proportional to |q|/M^, |p Ci d|/Mjv in the amplitude. Spin observables in pp reaction are 
universal and do not depend on the specific model. So, for simplicity we first consider the 
OBE amplitude for the process la. 

In the threshold limit the meson propagators and form factors become constants because 
they depend on the same variable k 2 — > k% = —M N M^ and using the explicit form of the 
functions II, E and W (cf. Appendix) we get 



n mdmft = \ .Uy.U.i. (2w 6 ) cos 9 b 6 mdmb , 

^m c m a = (! - ^ P )\jM N M^ (2m a ) cos 6 a 6- mamc , 

S m c m a = *(1 - «p) V M nM$ COS 6 a 5^ manlc 

W r mcma = ~(1 - k p )MI^M n {AM n + M„) {el! + (2m a )te; r } 5. maUlc . (24) 
TP' Z is not operative here, and one can express the amplitude in the simple form 



9 



with 



K(ko 



T T [ab] cd] = T U r [ab; cd] 



T = (1 - k p ) K{kl) M N Mf 2 {M N + ^M^ 2 , 
2 x 9nNN QpNN 9^ A 2 - m 2 A 2 - ml Af_ - ml Alt - ml 



P <Pf™ P ' 



{kl - ml){kl - ml) A 2 - k 2 Al - fcg AjL - k 2 A;% - P ' 



U r [ab;cd] = -A\/2m a m h cos6 h 8„ mam J mbmd d\ mar {6). (25) 

The spin is transferred to the meson only in the NNp vertex by a nucleon spin-flip, and 
for T ±1 {9 = 0) 7^ only for m a = ±~. Summation of |T r | 2 over the spin projections m c ^ of 
the outgoing particles leads to 

H I^qI 2 = 2T ^m a m b \ d 2m a r(9)\ 2 i ( 26 ) 

where A^ amt = 8 5 manib , and &Z a m b = 2 (1 + 8 5 maTO J. Using these equations we get the 
threshold limits for the beam - target asymmetry Cf^ = 1 for pp collisions. The diagonal 
spin density matrix elements are p^ 1 = | sin 2 9 and p±i±i = hO- cos2 ^)- Note that, as 
we mentioned above for the pp reaction, this predictions are universal because they reflect 
spin and parity conservation for two identical fermions and do not depend on the specific 
mechanism. All of the above discussed amplitudes are satisfy the same relations. 

For the pn reaction this kind of predictions depend on the model and generally on model 
parameters, say, on the relative contribution of neutral and charged meson exchange to 
the two-body T matrix in <p bremsstrahlung, shake-off and others. However, for the OBE 
amplitude la the threshold predictions for the beam-target asymmetry and spin-density 
matrix are parameter independent. From Eq. (|26|) we get for the beam-target asymmetry 
= 0.8, which means that the ratio of singlet to triplet spin states is ~. The diagonal 
spin density matrix elements are the same as in the pp reaction. 

Using Eq. we can make an important prediction: the ratio of the unpolarized cross 
sections in pn and pp reactions is a np /a PP = 5/2 ■ \ = 5, where the factor \ in the pp cross 
section represents the symmetry factor. This prediction does not depend on the OBE model 
parameters. As we shall see below, the OBE amplitude la dominates, hence fore we can 
expect that the ratio of the corresponding total cross sections would be close to the factor 
5. 



VI. RESULTS 

Our results for the total cross section of the pp reaction are shown in Fig. 2 by the solid 
lines. As we mentioned above, the comparison of the momenta at the 7rp0 vertex in the 
AW — > NNcj) reaction with the irN — > Nip reaction shows a limited reliability of the pure 
OBE model to the region of A s = yfs — 2M^ — < 0.15 GeV. One can easily extend 
the model to a wider energy region, say A s ~ 1 GeV, by reducing the cross section with 



10 



the effective cut-off factor oc exp(— 0.8A S ) resulting from the energy depending coupling 



strengths in Ref. |T(|. The resulting cross section is shown by the dashed lines. 

Triangles in Fig. 2 show our prediction for the <fi bremsstrahlung contribution lb. Open 
and full circles show the predictions for the shake-off channel lc with a strangeness proba- 
bilities B 2 = 0.01 and 0.05, respectively. The dot-dashed lines show the contributions from 
the pure OBE channel la. The <fi bremsstrahlung and shake-off contributions have the same 
order of magnitude. The shake-off channel in this reaction is reduced strongly by the form 
factor V in Eq. ( p2]) . The solid lines show the coherent sum of all amplitudes with B 2 = 0.05. 
The left (right) panel in Fig. 2 displays the result for a purely real (imaginary) phase of the 
two-body T matrix. The interference between the pure OBE amplitude la and the shake-off 
amplitude lc in the first case is rather small and therefore the difference between the results 
of the shake off with a\ = ±1 is not seen. But it may be seen, in principle, in second 
case, where the positive value of a\ is marked by crosses. In spite of the difference between 
the separate contributions of the conventional OBE channel la and the exotic <fi production 
channels, which is more than an order of magnitude at A s ~ 0.1 GeV for the above given 
parameters and g^NN — —0.24 and B 2 < 0.05, the interference may amount up to 50% of 
the total cross section. 

The cross section for the pn reaction is greater by a factor ~ 5 as given already in our 
threshold-near prediction. 

The interference between the conventional OBE process la and the other exotic <ft chan- 
nels lb,c is much stronger in polarization observables. Fig. 3 shows our finite-energy predic- 
tions for the beam - target asymmetry Eq. (1) as a function of the (ft production angle 9 at 
fixed recoil nucleon angles 9 p ',ip p > = 7r at A s = 0.1 GeV and |q| = §A(s, M|, 4M^ r )/(2y / s) 
for the pp reaction for separate channels (A is the usual triangle function). The contribution 
from the OBE channel la is depicted by the dot-dashed line, while the bremsstrahlung 
contribution lb is marked by triangles, and the <j) shake-off process lc is marked by dots. 
The threshold prediction is the dashed line. 

One can see that the asymmetry for the conventional OBE channel la and the <fi 
bremsstrahlung channel lb are qualitatively very similar. But in the case of a real phase of 
the two-body T matrix there is a strong "destructive" interference in the asymmetry which 
forces a significant decrease of the total asymmetry. This is shown in Fig. 4 (left panel) where 
we display the beam - target asymmetry as a coherent sum of all channels. Solid, dotted 
and dash-dotted lines correspond to B 2 = 0., 0.01 and 0.05, respectively. The contribution 
from the shake-off process is rather small. The strong decrease of the total asymmetry may 
be considered even as interesting quantitative manifestation of the bremsstrahlung pro- 
cess lb. If one would use instead of q^nn = —0.24 the same positive quantity, then the 
asymmetry results in much larger values 0.4 ■ • • 0.6; the corresponding result is marked by 
crosses. 

The right panel of Fig. 4 illustrates the asymmetry for a purely imaginary two-body T 
matrix. The crosses mark a positive value of a x in Eq. ([19|) . The sign of a± is still unknown. 
The area between marked and unmarked dot-dashed lines represents the theoretical uncer- 
tainty of the predictions for B 2 = 0...0.05. For positive values of ai the asymmetry differs 
noticeably from the pure OBE process la prediction and may reach even negative values. 

The influence of this uncertainty is smaller for the spin density matrix Eq. (^]). Fig. 5 
(left panel, same notation as in Fig. 3) depicts p o separately for each channel at A s = 0.1 
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GeV. The anisotropy of the angular distribution in 2-body decays — > a + a~ is given by 
Eqs. (0,0). One can see in Fig. 5 (right panel, same notation as in right panel of Fig. 4) 
a strong modification of the anisotropy for — > e + e~ decays by the shake-off mechanism at 
small angles 9. This effect is smaller for K + K~ decays. The dotted and dot-dashed lines 
marked by diamonds correspond to calculations with real two-body T matrix and B 2 =0.01 
and 0.05, respectively. 

In pn reactions all these polarization effects exist too, but they are not so strong as in 
the pp reactions. That is because the OBE channel la here is a factor of 5 greater, while 
the contribution of the "exotic" channels remains on the same level as in the pp case. 

Finally, we would like to mention that the choice A^ p7r = 2.9 GeV (instead of 0.77 GeV) 
increases the OBE contribution la by factor 4 • • • 6, while the other "exotic" contributions 
remain on the same level. This results in a decrease of their relative role in polarization 
observables which then become much closer to the pure OBE channel la predictions. 



VII. SUMMARY 



In summary we calculate within an extended OBE model with "exotic" production 
channels the cross section and polarization observables for the reaction NN — *> NN(f). We 
predict the threshold behavior of cross sections for polarized and unpolarized reactions. 
Above the threshold the polarization observables differ strongly from the threshold pre- 
dictions and behave differently for each channel. This means that measurements of the 
beam - target asymmetry and anisotropy of the decays may reveal the presence of hidden 
strangeness and help disentangling various reaction mechanisms. 
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members of the DISTO collaboration are gratefully acknowledged. One of the authors 
(A.I.T.) acknowledges the warm hospitality of the nuclear theory group in the Research Cen- 
ter Rossendorf. This work is supported by the BMBF grant 06DR829/1 and the Heisenberg- 
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APPENDIX A: FUNCTIONS II, S AND W FOR THE OBE AMPLITUDE 

^-m d m b (Pb,Pd)- 

Using the definition of the Dirac spinor 



and the identity 

Xf (cr ■ A) Xi = 2miA z 5 mimf + (A x + 2imiA y )6 mi - mf , (A2) 

we find 

^-m d m b = C bd [ 2m d (a b cos 9 b - a d cos 9 d ) 5 mbmd 

- a d sin 9 d exp(-2im d y?d) 5 m6 _ md ] , (A3) 
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where Cm = ^{E b + M N )(E d + M N ). When the incoming particle moves along the z axis, 
cos 9 h = ±1. 



^m e m a (^p) : 

Calculating Yl u (k p ), where k p = p c — p a , we use Gordon's identity 

u(Pc) [Y] U(p a ) = ^HPc) [{Pa+PcT + 1(7^ k Pl/t ] u{p a ). 

which allows to express T, u as the sum 

yu yu I yu 



yu 

^2 m c m a 



(1 - Kp)M(p c )7 !y M(p a ) 



with P ac = p a + p c - Direct evaluation of Si gives 



yu 



2M; 



-P^ c C ac [(l - a a a c cos 6> a cos 6> c ) 5 manic 



N 



2m a a a a c cos 9 a sin C exp( — 2ifn c ip c ) S ma ^ r 



Explicit expressions for the components of S2 read 



J 2ac 



y y 



J 2ac 



K p )C ac [ (1 + a a a c cos 9 a cos 9 c ) 5 mamc 

+ 2m a a a a c cos 9 a sin 9 C exp(—2im c (p c ) 5, 
K p )C ac [ a c sin 9 C exp(-2im c ip c ) 5 mamc 

+ 2m a (a a cos9 a - a c cos9 c ) <5 ma _ m J, 
% (1 - k p ) C ac [ 2m a a c sin9 c exp(2im c Lp c ) 5 mainc 

+ (a a cos 9 a - a c cos 9 C ) 8 ma - mc ], 
-- {I - Kp) C ac [ (a a cos 9 a + a c cos 9 C ) 5 marHc 

+ 2m a a c sin9 c exp(-2im c ip c ) 5 mo _ m J. 



m a -m c \i 



W^ c (p a ,p c ): 

For calculating W r we use the identity 

e^aptf = k pG q • (S x - q k p • (S x s* r ) 

+ S e* r • (k p X qr) - £* r S • (fc p X qr), 

which allows to express amplitude W r in compact form as 

W r —i (A - C r + B ■ D r ) , 

where 

A = k p0 q - q k p , B = k p x q, 
C r = £ x e* r D r = S e* r - e* r E. 



(A4) 



(A5) 



(A6) 



(A7) 
(A8) 
(A9) 
(A10) 



(All) 



(A12) 



(A13) 
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(for simplicity we skip here the arguments p a , p c , ... in W r and S r ). The independent 
meson polarization vectors e r = (e r x , e r y1 e r z , Eq) with r = 1, 2, 3, read 

e 1 = (cos6, - sin#, 0), 
e 2 = (0, 1, 0, 0), 

£ 3 = (7 </ ,sin6», 0, 7 cos6», 7^), (A14) 
where 7^ = E^/M^, = \q\/E^,. Using the polarization vectors in the helicity basis 

e *° = £ 3, e^^T^Tie 2 ) (A15) 

we get the amplitude in the helicity basis as 

W° = W 3 , W ±l = T-^iW 1 TiW 2 ), (A16) 

v2 

which may be directly used in numerical calculations. 
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FIG. 1. Diagrammatic representation of processes of the (j) production: emission from (a) the 
internal meson line, (b) the external nucleon legs, and (c) the exchanged meson - nucleon vertex. 
Exchange diagrams are not displayed. 
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FIG. 2. The total cross section of the pp — > pp0 reaction as a function of the excess energy. 
Data (fat dots) from j27j. The individual contributions are shown separately: pure OBE process 
la (dot-dashed lines), <j) bremsstrahlung lb (triangles), 4> shake-off (open and filled circles for B 2 = 
1 and 5%). Solids lines depict the coherent sum of all amplitudes (with i? 2 =0.05), while dashed 
lines represent these sums with additional energy depending effective cut-off described in text. Left 
panel: real two-body T (the difference between positive and negative values a\ is not seen on this 
scale); right panel: purely imaginary two-body T matrix (crosses mark the coherent sum with 
positive value of a\). 
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FIG. 3. The finite - energy beam-target asymmetry for the pp — > ppcfi reaction as a function of 
the c.m.s. <j) production angle 9 for different channels. The contribution from the OBE channel la 
is the dot-dashed line, the (f> bremsstrahlung channel lb is marked by triangles, and the <p shake-off 
lc is marked by dots. The threshold prediction is depicted by the long-dashed line. 
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FIG. 4. The beam-target asymmetry as a coherent sum of all channels. Solid, dotted and 
dot-dashed lines correspond to B 2 =0., 0.01 and 0.05, respectively. Left panel: real two-body T 
matrix (the difference between positive and negative a\ values is not visible on this scale; crosses 
indicate the result for g^NN = +0.24); right panel: imaginary two-body T matrix (a positive value 
of cti is marked by crosses). The threshold prediction is depicted by the long-dashed line. 
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FIG. 5. The spin density matrix element poo separately for different channels (left panel; 
notation as in Fig. |3|) and the asymmetry of (f) — > e + e~ decays normalized to the threshold value 
(right panel; for an imaginary two body T matrix; notation as in right panel of Fig. |I|; the results 
for a real phase are marked by diamonds). 
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